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ABSTRACT 

We employ a low angular momentum accretion-outflow scenario to model the flares 
emanating out from the central region of Sgr A*. The primary donor for matter ac- 
creting onto the central SMBH of Sgr A* is assumed to be the WR star ISR 13 E3. 
We analytically calculate the specific energy and angular momentum density of stel- 
lar wind originating from ISR 13 E3 and study the dynamics of that wind down to 
the very close vicinity of the central SMBH of Sgr A*. We show that on the way to 
the Galactic centre, such wind-fed accretion may encounter standing shocks and such 
shock drives outflow from the close vicinity of the SMBH. Matter content of such out- 
flow is computed and it is argued that such outflow is responsible for production of the 
Galactic centre flares. We then self-consistently compute the luminosity Cj and the 
duration time scale t of such flares, as a function of fundamental accretion parame- 
ters. Our theoretical calculation of Cj and r are in good agreement with observational 
results. 
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1 INTRODUCTION 

The center of our Galaxy harbors a massive black hole, and 
the surrounding region including the central SMBH is now 
customarily referred as Sgr A* as a whole, after the radio 
source first discovered at that location (for a review, see e.g. 
Melia & Falcke 2001). Sgr A* shows frequent flares originat- 
ing from the direct vicinity of the central black hole. In the 
X-ray band, multiple fiares are superimposed on a steady, 
extended emission at the level of ~ 2.2 x 10'^'^ erg s~^ cm~^ 
(Baganoff et al. 2003a). Two extremely bright flares have 
been reported so far (Baganoff et al., 2001, Porquet et al. 
2003; maximum flux of 1.0 ± 0.1 x 10^^ and ?>.&t°ol x 10^^ 
respectively), along with many fainter flares observed in the 
Chandra data (e.g. Eckart et al. 2004). The duration of the 
flares ranges from half an hour to several hours, while the 
rise/decay time is found to be of the order of few hundred 
seconds (Baganoff 2003). 

Variable emission was also detected in the NIR band. 
Quiescence emission takes place at the level of ~ 1.9 mjy 
(Eckart et al. 2004). Variable and quiescent emission was 
reported by Genzel et al. (2003) based on the VLT obser- 
vations, and by Ghez et al. (2004) based on Keck data. In 
two of the events, a 17 min periodicity was found (Genzel et 
al. 2003). X-ray and NIR outbursts are directly related, as 
shown by the detection of a simultaneous NIR/X-ray event 



(Eckart et al. 2004). The 2-8 keV luminosity of the event was 
~ 6 X 10^^ erg s"^ cm"^ and 2.2 /im flux was 3.7 mJy. The 
duration of the X-ray event was 55 - 115 min, with the de- 
cay time of such flares measured to be of the order of several 
minutes. Detection of TeV emission from the Galactic Cen- 
ter (Aharonian et al. 2004) supports the view that emission 
comes from a jet. 

These properties of the flare emission suggest that such 
flares originate in the innermost region of accretion flow 
onto the central black hole, and we need to introduce an ap- 
propriate theoretical model of accretion to understand the 
properties of these flares. Most of the works on accretion 
processes onto the SMBH of Sgr A* concentrate either on 
almost purely Bondi (1952) type accretion, or on high angu- 
lar momentum ADAF type flow (Narayan, Yi & Mahadevan 
1995; Yuan, Quataert & Narayan 2004), which is sometimes 
claimed to be coupled with jet-like outfiows (e.g. Markoff et 
al. 2001, Yuan, Markoff & Falcke 2002). 

In this letter, however, we employ a different kind of 
accretion model - a low angular momentum (highly sub- 
Keplerian) flow with standing shock, to explain the gener- 
ation of flares from Sgr A* and to compute the luminosity 
of the flare Cj along with its duration time scale r in a 
self-consistent way. Sub-Keplerian advective accretion onto 
galactic and extra-galactic black holes is expected to pro- 
duce multi-transonic behaviour (see e.g. Das 2004, Barai, 
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Fig. 1: The schematic picture of the wind flow from the star IRS 
13E3 to-wards the Galactic center. 

Wiita & Das 2004, and references therein) and standing 
shocks are an essential ingredient in an low angular mo- 
mentum flow in general (Das 2002, hereafter D02, Das, 
Pendharkar & Mitra, 2003, hereafter DPM, and the refer- 
ences therein). Such shocks in turn play an important role 
in governing the overall dynamical and radiative processes 
taking place in accreting material and the hot, dense, exo- 
entropic post-shock fluid is supposed to be responsible for 
launching jets/outflows from black hole accretion discs (Das 
& Chakrabarti 1999, hereafter DC, Das, Rao & Vadawale 
2003, hereafter DRV, and references therein). Such a cou- 
pled accretion-outflow model may reproduce the character- 
istic behaviours of the Sgr A* flares in a natural way. 



2 SOURCE OF THE ACCRETING MATERIAL 

SteUar winds originating from the central cluster are plen- 
tiful sources of gas. Most of this material is likely to be 
expelled from the central region (e.g. Quataert 2003) but a 
remaining fraction may power the observed activity. Among 
various mass loosing stars, particularly active one is the 
Wolf-Rayet star IRS 13 E3 (Paumard et al. 2001, Melia & 
Falcke 2001). We assume this star to be the dominant source 
of the matter accreting onto Sgr A* . We determine the angu- 
lar momentum and the Bernoulli constant of the flow from 
IRS 13 E3 in the following way: 

We consider the case where the wind velocity is much 
higher than the orbital velocity of a star. Therefore, the 
material ejected from a fractional region of the star surface 
located at 4>o (see Fig.l) can reach the gravity center with 
zero angular momentum (see e.g. Loeb 2004). This angle is 
given by the condition: 

, ^star / ^\ 

Sm0o = . (1) 

^wind 

The wind is not signiflcantly accelerated as long as the 
gravity field is dominated by the stellar component. It is 
moderately supersonic, with the Bondi-Hoyle-Lyttleton ac- 
cretion radius, -Rbhl, given by the formula 

2GM 



(2) 



where tiwind is the flow velocity and Vs is the sound speed. We 
further assume that the velocity t^wind is unmodified down 
to a distance -Rbhl. The value of the Bernoulli constant, £, 



£ = 0.5t;^ind H ^ 

7i - 



GM 



1 ^Jbhl 



(3) 



is effectively determined by the value of the polytropic in- 
dex on the inflow, 7^, and an assumed gas temperature. 



where rg = 2GAfBH/c^. Through out this work, we use the 
Paczyiiski and Wiita (1980) potential to describe the flow. 

The spherically-symmetric wind blowing a.t cj) ^ (j>o or 
out of the orbital plane will posses certain amount of positive 
or negative specific angular momentum (angular momentum 
density), A. In the second order approximation 



A ^ — Uwind-Di5<^cos(<^o)[l — l/2tan(j!)o5<j!>]. 



(4) 



where 50 = — ^o, (j) is the azimuthal angle of the element 
at the star surface and D is the distance between the star 
and the Galactic center. 

A cylindrical fraction of this fiow, with = Rbhl/D, 
and A6 ^ -Rbhl / D will be intercepted by the central black 
hole, where the angle 9 determines the deviation from the 
orbital plane. 

Integrating the Eq. 4 with respect to S(j} and 56 in the 
limits specified by A(j) and A9, we obtain the net angular 
momentum of the flow as: 



Ac 



/ Rbhl \ ' 



(5) 



\/l - (Ustar/fwind)^ 

The above relation is valid if the star velocity is significantly 
smaller than the wind velocity. Large value of the angular 
momentum density of the donor star, v^ti^D, is decreased 
by small quadratic term in the -Rbhl/I? ratio. 

The formula changes if the wind is not perfectly 
isotropic. If the departure of the mass flux, 5(pwind«wind)/ < 
Pwindi'wind >, happen at distances (5<^-Rstar, the net angular 
momentum density can be estimated by integrating the for- 
mula 4 as 



Ac 



1 / -RbHl\ <5(Pwind1'wind) 

--DUstar r, — -, 

TT \ U J < Pwindfwind > 



(6) 



i.e. the net angular momentum density is a fraction of 
i?BHLWwind- The flow, initially one-sided, becomes roughly 
spherical below -Rbhl- We further express A in dimension- 
less units 2GM/c, and the radius in rg. 



3 LOW ANGULAR MOMENTUM FLOW 
CLOSE TO SGR A* 

If the accreting material is assumed to be at rest far from 
black hole, the fiow must exhibit transonic behaviour in or- 
der to satisfy the inner boundary conditions imposed by 
the event horizon. Low angular momentum fiow may posses 
more than one sonic point, as first shown by Abramowicz & 
Zurek (1981). Typically the external sonic point, rout, lies 
close to the corresponding Bondi radius. The internal sonic 
point rin and the middle sonic point rmid exist within and 
outside the marginally stable orbit, respectively, for general 
relativistic (Das 2004, Barai, Das & Wiita 2004) as well as 
for post-Newtonian (D02, DPM) model of accretion flow. 
The location of the sonic points can be calculated as a func- 
tion of the specific flow energy £ (the Bernoulli's constant), 
angular momentum A and inflow polytropic index 7^ (see, 
e.g. §3 of D02 for details of such calculations). 

If A is almost zero, a shock does not form, and accre- 
tion remains supersonic down to the event horizon after it 
crosses rout. For slightly larger A the centrifugal barrier be- 
comes strong enough, infiowing matter starts pilling up close 
to the black hole due to the resistance offered by the bar- 
rier, and the depleted matter may break the incoming fiow 



© 0000 RAS, MNRAS 000, 000-000 



Low angular momentum accretion- outflow model of flares from Sgr A* 3 



behind it and consequently a shock forms. Such shocks may 
become steady and standing so that they can be studied 
within the framework of stationary flow (D02, DPM and 
references therein). 

Following D02, wc consider here a stationary, non-self- 
gravitating, non-magnetized, inviscid accretion of polytropic 
fluid. We assume that the flow proceeds through a standing 
shock, so the discontinuity in the radial velocity allows to 
match the supersonic flow through rout with the subsonic 
flow through rin. The exact location of the shock as a func- 
tion of parameters [f,A,7f], is obtained by solving the gen- 
eralized Rankine-Hugoniot conditions. We assume that the 
shock is nou-radiating and infinitcsimally thin. 

Formation of the shock leads additionally to the gen- 
eration of an outflow, and the shock location, r^h, may be 
assumed to be the radial length scale of the outflow launch- 
ing zone. The exact amount of the outflowing material Mout 
can be calculated as a function of [£", A, 7^, 70], where 70 is 
the outflow polytropic index and is always less that 7i due 
to the radiation momentum deposition on slowly expanding 
outflowing matter at shock location. In our work, we as- 
sume 7o to be a free parameter (subjected to the constraint 
7o < 7i)) although in reality 70 may be directly related to 
the heating and cooling processes taking place in post-shock 
matter; see DC and DRV for details. 

Flares observed from Sgr A* can be due to the emission 
from this outflowing material. Jet emission seems to be the 
plausible origin of both the X-ray and the NIR radiation 
(e.g. Markoff et al. 2001), if the dissipation within the shock 
itself does not lead to a strong emission. In this case the 
luminosity of Sgr A* during the flare, £j, is related to the 
outflow rate, Mout, as Cj = O.lc^Afout. Here we assume the 
10% efficiency of energy conversion, since studies of AGN 
jet indicate rather high efficiency of jets even if the radia- 
tive efficiency of the accretion flow is low (e.g. Maraschi & 
Tavccchio 2003). Wc basically model a stationary situation. 
However, the flow pattern with an outflow can fully develop 
under the condition that the accretion phase with a shock 
lasts long enough. Therefore, the minimum duration of the 
flare, r, is associated with the infall time scale of post-shock 
accretion and can be defined as: 



dr 
u{r) 



(7) 
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Fig. 2: Example of a typical flow topology with a shock. The 
position of the shock, r^h. , is determined by the Rankine-Hugoniot 
conditions. 
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where u{r) is the dynamical flow velocity and is the loca- 
tion of the event horizon. 

It is to be noted that multi-transonic flow and standing 
shocks form for a specific region of parameter space spanned 
by [£, A, 7i] (sec, e.g. Figure 4. of D02 for a global classifica- 
tion of shock forming parameter space) . For certain values of 
[£, A,7j], no multi-transonic or standard mono-transonic sta- 
tionary solution exists. The corresponding flow pattern be- 
comes inherently time-dependent. Such non-stationary flow 
might even better account for Sgr A* variability but we have 
no adequate description of such a flow. 



4 RESULTS 

We adopt 3.8 x IO^Mq for the black hole mass in Sgr A* 
(Ghez et. al. 2003) and consider the accretion rate to be 
equal to 10~*Mo yr~^. The distance to the donor star is 
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Fig. 3: The relation between the shock location and the outflow 
rate. Marked along the curve are the corresponding values of the 
angular momentum density in 2GM/c units. 

3.5 pc, and the wind velocity of the star, Vwind, is of order 
of 1000 km/s (adopted after RockefeUer et al. 2004). We 

assume that the orbital star velocity, Vstai, is 200 km/s, of 
order of radial velocities of other stars at similar distance 
from the center although the actual radial velocity of this 
star is lower (Eckart & Genzel 1997, Paumard et al. 2001). 
As the polytropic index of the inflow, we take the value 71 = 
1.426, representative for non-relativistic flow, and for the 
wind temperature we take 1.3 kcV, the plasma temperature 
estimated at the basis of the analysis of the extended X-ray 
emission (Baganoff et al. 2003a). This means that the flow 
is moderately supersonic at the 7?bhl (vs = 600 km s~^, 
-Rbhl = 7.4 X 10^^ cm or 6.6 x lO'^rg). These values allow us 
to estimate the Bernoulli constant of the flow, £, in Eq. 3 as 
5.2 X 10~® in units of c^. The effective angular momentum 
given by Eq. 5 is usually small, Aes = 0.14. This means that 
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typically the flow can proceed directly toward the black hole, 
without the need for angular momentum loss and without 
considerable dissipation. However, if the wind is occasionally 
non-uniform, 5(pwind'ywind)/ < Pwind"wind > of order of 2 per 
cent, changes the not angular momentum by a factor of 11. 
Such a large value of the angular momentum forms effective 
centrifugal barrier and the inflow may proceed through a 
shock, with accompanying outflow and a burst of radiation 
in the X-ray and NIR bands. 

Fig. 2 represents a characteristic topology of shocked 
accretion flow. The values of [f,A,7i,7o] used are provided 
in the figure. Matter first passes through rout ( 3.2 x 10*rg) 
and encounters a shock at r^h (10.03 r-g) close to the black 
hole. The dashed vertical line marked with a down-ward ar- 
row represents the shock transition. M_ and M+ are the 
pre/post shock Mach numbers and the shock strength <S is 
defined as 5 = j^-, which comes out to be 16.64 for this 
case. Post-shock subsonic inflow becomes supersonic again 
after crossing rin (2.61 rg) and finally dives through the 
event horizon T e. Part of the shock compressed hot and dense 
matter with polytropic index 70 (=1.327) emerges as out- 
flow. The outflow rate Mout in this case comes out to be 
2.1 X 1O~^°M0 yr~^. The corresponding values of Cj and 
r comes out to be 1.1x10^® erg s~^ cm~^ and 4 x 10^ s 
respectively. 

The position of the shock is quite sensitive to the choice 

of A, if other parameters of the flow are kept fixed. One 
can have a range of r^h and Mout by varying [£,A,7i,7o] 
as well as by varying A only by keeping [f,7i,7o] constant. 
For a fixed set of ,7i, 7o] shown in the figure, we represent 
the dependence of the outflow rate on the location of the 
outflow launching zone in Fig. 3 by varying the value of flow 
angular momentum. Similar figure can be drawn for other 
appropriate values of [f ,7i,7o] as well. Note, however, that 
although very large values for rgh (outflow launching zone) 
can be obtained as a consistent mathematical solution, they 
may not correspond to the real physical situation. Smaller 
values of the outflow launching zone length scale could be 
obtained by further decreasing the value of 7i. 

Fig. 4 shows the relation between the luminosity of the 
flare Cj (in CGS unit) and the flare duration time scale r (in 
units of one thousand seconds) . The figure is drawn for the 
same [f ,7i,7o] used to draw the Fig. 3. We show the values 
of Cj and t for the range of r^h as r^h = 8.69 Vg to r^h = 17.03 
Tg, the corresponding range for A is from 1.651 to 1.675. 
The low luminosity flares axe of shorter duration and the 
length scale at which they are generated is also shorter, this 
is because the dynamical time is shorter as matter gets closer 
to the Black hole. One can have a further small (compared 
to the lowest value shown in the figure) value of r by fine 
tuning the values of [A,7i,7o] for a fixed value of £. Note, 
however, that r co-relates with the black hole mass. If one 
considers Mbh = 2.6 x 10*^ Af© (the lower limit of Mbh 
for Sgr A*), all values of r will be reduced. If T2.6 and ra.g 
corresponds to Mbh = 2.6 x lO'^M© and Mbh = 3.8 x 
W^Mq respectively, then (rs.s — T2.6) comes out to be about 
1250 seconds for low luminosity flares and the difference 
increases for flares with higher value of luminosity, i.e., for 
larger value of Vsh- For example, for Vsh = 17.03rg, T2.6 is 
only equal to ~ 7449 seconds whereas T3.8 comes out to be 
equal to ~ 10887 seconds. 
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Fig. 4: The relation between the flare luminosity Cj and flare 
duration time scale r (in thousands of seconds) for some repre- 
sentative cases, see text for detail. 

5 DISCUSSION AND CONCLUSIONS 

The type of weakly rotating flows presented in this letter, 
have not been theoretically explored much in the literature; 
although thoy arc well exhibited in nature for various real 
physical situations like detached binary systems fed by ac- 
cretion from OB stellar winds (lUarionov & Sunyaev 1975; 
Liang & Nolan 1984), semi-detached low- mass non-magnetic 
binaries (Bisikalo et al. 1998) and super-massive BHs fed by 
accretion from slowly rotating central stellar clusters (II- 
larionov 1988; Ho 1999 and references therein). Even for 
a standard Keplerian accretion disc, turbulence may pro- 
duce such low angular momentum flow (e.g. Igumenshchev 

6 Abramowicz 1999, and references therein). 

We consider this type of flow as an attractive model 
of Sgr A* activity. We expect that accretion onto the cen- 
tral black hole proceeds roughly continuously but small (up 
to 2-3 %) fluctuations in the wind density or velocity lead 
to small variations in the angular momentum density which 
in turn results in temporary shock/jet formation and con- 
sequently an enhanced dissipation. The range of predicted 
shock positions lOrg < rah < 300rg give interesting range 
of expected burst timescales, and more luminous flares are 
expected to last rather longer than fainter flares. 

The predicted exemplary flare luminosities are in the 
interesting range. The exact estimate of the flare bolometric 
luminosity is rather difficult, as nicely discussed by Aharo- 
uiau & Ncronov (2004). As on order of magnitude estimate 
we can assume that the bolometric luminosity of the flare is 
ten times higher than the measured X-ray flux so the pre- 
dicted outburst would have the measured X-ray luminosity 
around 10'^" erg s^^ cm^^. 

The shown luminosity range is, however, rather narrow. 
The barionic load of the outflow (and hence, Cj) is controlled 
mainly by the post shock thermal pressure and post shock 
thermal energy generated is modulated by the total energy 
content of the flow as well as the thermal energy content of 
fluid. Hence broader luminosity range can be obtained if we 
allow for variation of 7i,7o and £ on the top of perturba- 
tions of the angular momentum density. We found that Cj 
correlates with [£,7^,70]. While the variation of Cj is quite 
sensitive on 70, it is relatively less sensitive to the variation 
of [£, 7j] . Hence the thermal properties of the outflow mainly 
contributes to the variation of Cj as is expected. 
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Our approach to description of the low angular momen- 
tum flow is based on several assumptions. 

First, in our work, viscous transport of the angular mo- 
mentum is not explicitly taken into account. Even thirty 
years after the discovery of standard accretion disc the- 
ory (Shakura & Sunyaev 1973), exact modeling of viscous 
multi-transonic BH accretion, including proper heating and 
cooling mechanisms is still quite am arduous task. Neverthe- 
less, extremely large radial velocity close to the BH implies 
Tinf << Tiiisc {Tinf and T^,isc arc the fall and the viscous time 
scales respectively). Large radial velocities even at larger 
distances are due to the fact that the angular momentum 
content of the accreting fluid is relatively low (sec, e.g., Be- 
loborodov & lUarionov 1991, Igumenshchev & Beloborodov 
1997, Proga & Begelman 2003). Hence, our assumption of 
inviscid flow is not unjustified. 

On the other hand, the introduction of viscosity would 
further reduce the radial angular momentum of the accreting 
matter. As we have seen (Fig. 3), lower values of A produces 
the smaller values of the outflow launching zone (rsh) and 
smaller amount of Mout , which means (see Fig. 4 and related 
discussions) that viscous transonic flow would produce flares 
from Sgr A* with shorter duration time scale (r), compared 
to the values of r obtained using our inviscid flow treatment. 

Second, we apply here the stationary solution to deduce 
the properties of the time-dependent flow. However, explicit 
time-dependent considerations addressed so far only spe- 
cific issues like time-dependent and oscillatory behaviours 
of the shock (see, e.g. Okuda, Teresi, Toscano & Molteni 
2004, and references therein), which, for example, leads to 
the quasi-periodic oscillation of black hole candidates (Das 
2003, and references therein). Several groups performed nu- 
merical liydrodynamical computations. For example, numer- 
ical simulations by Coker & Melia (1997) and Rockefeller et 
al. (2004) of winds from many stars lead to predictions of 
much higher average angular momentum than adopted in 
our paper (A ^ 60) but also to much higher accretion rate 
than allowed by the Faraday rotation constraints (BaganofF 
et al. 2003a). The problem is that such computations can- 
not resolve the flow deeply inside Rbhl, where the flare 
formation occurs. 

Our picture of accretion flow is qualitatively similar to 
the one developed by Loeb (2004). However, he considered 
accretion from stars very close to the black hole, well within 
appropriate Bondi-Hoyle-Lyttleton radius, with less efficient 
winds and star velocities being smaller than the wind veloc- 
ity only at a fraction of an orbit. In this picture accretion 
from any of the considered stars is expected to occur once 
per orbit and to last only for a few months. IRS 13 E3 is 
expected to supply the mass continuously. Therefore, Loeb 
(2004) model would require sigruficant correlation of the long 
time-scale activity with the motion of nearby stars while our 
picture would be in agreement with no such trends. Obser- 
vational search for such correlations will help to determine 
the dominant source of the accreting material. 
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